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Abstract. This paper summarizes the information gathered for 16 still unpublished exoplanet candidates discov- 
ered with the CORALIE echelle spectrograph mounted on the Euler Swiss telescope at La Silla Observatory. 
Amongst these new candidates, 10 are typical extrasolar Jupiter-like planets on intermediate- or long- 
period (100<P< 1350d) and fairly eccentric (0.2 < e < 0.5) orbits (HD 19994, HD 65216, HD 92788, HD 111232, 
HD 114386, HD 142415, HD 147513, HD 196050, HD 216437, HD 216770). Two of these stars are in binary systems. 
The next 3 candidates are shorter-period planets (HD6434, HD 121504) with lower eccentricities among which 
we find a hot Jupiter (HD 83443). More interesting cases are given by the multiple-planet systems HD 82943 and 
HD 169830. The former is a resonant P2/P1 =2/1 system in which planet-planet interactions are influencing the 
system evolution. The latter is more hierarchically structured. 
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1. Introduction 

. For more than 5 years the CORALIE planet-search pro- 
' gramme in the southern hemisphere llUdrv et al.ll2nn0a^ 
] has been ongoing at the 1.2-m Euler Swiss telescope - 

■ designed, built and operated by the Geneva Observatory 
] - at La Silla Observatory (ESO/Chile). During these 5 

■ years, CORALIE has allowed the detection (or has con- 
] tributed to the detection) of 38 extra-solar planet candi- 
' dates. This substantial contribution together with discov- 
eries from various other programmes have provided a sam- 
ple of more than 115 exoplanets that now permits us to 
point out interesting statistical constraints for the planet 
formation a nd evolution scenarios (see e.g. iMavo mom 
lUdrv et alJl2003d: .Santos ct al. 2003b- for reviews on dif- 
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* Based on observations collected with the CORALIE echelle 
spectrograph on the 1.2-m Euler Swiss telescope at La Silla 
Observatory, ESO Chile 

** The precise radial velocities presented in this paper 
are available in electronic form at the CDS via anony- 
mous ftp to cdseirc.u-strasbg.fr (130.79.128.5) or via 
|http : // cdsweb . u-strasbg . f r/ cgi-bin/ qcat ? J/ A+A/ except 
for the multi-planet systems that will appear in a future pa- 
per describing their dynamical evolutions, taking planet-planet 
interaction into account 



ferent aspects of the orbital-element distributions or pri- 
mary star properties). 

The majority of our CORALIE exoplanet candidates 
have been published in a series of dedicated papers^, the 
latest among them reporting the detection of the shortest- 
period Hot Jupite r discovered by ra dial-velocity surveys 
around HD 73256 llUdrv et al.ll2003c|l and the very inter- 
esting case of HD 10647 ljUdrv et al.ll2003a|) . a star with a 
high IR excess indicative of the presence of a debris disk. 
The present paper of this series describes the CORALIE 
exoplanets that have not been published yet. This subsam- 
ple includes candidates announced several months ago, 
rapidly after their detection to allow follow-up observa- 
tions. It also includes some candidates with very long pe- 
riods or that arc members of multi-planet systems requir- 
ing a delay in their final analysis. Also, some of the new 
candidates correspond to very recent detections. 

The paper is organized as follows. In the next section 
we summarize the primary star properties. The radial- 
velocity measurements and inferred orbital solutions will 



^ Another dedicated series has also been started for close bi- 
naries requiring a 2-dimensional correlation analysis for radial- 
velocity estimates. The method ha s already revealed a 2.5- 
Mjup planet orbiting the primary JZucker et alJ l2003b^ and 
a 19-Mjup brown dwarf orbiting the secondary 1 Zucker et alJ 
1200331" ) of the HD 41004 close visual binary system 



2 



M. Mayor et al.: The CORALIE survey for southern extra-solar planets. XII 



Table 1. Observed and inferred stellar parameters for the stars hosting planets presented in this paper. Definitions 
and sources of the quoted value s are given in the text. The age an d rotational period estimates are based on calibrations 
of the R'uj^ act ivity indicato r llDonahu e 1993: Noves et al .''l984V whose ref erence source is also indicated: (S) for this 
paper foUowing lSantos et alJ l|2000() . CH) for .Henrv et all H199a and (B) for lButler et~al] l)2002l) . The apphed analyses 
and uncertainty estimates can be found in the quoted references. 
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be presented in Sect. 3. In the last section we summarize 
the results and provide some concluding remarks. 



2. Parent star characteristics 

The CORALIE planet-search t argets have been selected 
from the Hipparcos catalogue (^ES3[1993). Our planet- 
star subsample benefits thus from the photometric and 
astrometric information gathered by the satellite. 

A high-resolution spectroscopic abundance study was 
performed for most of these stars by N.C. Santos in 
his study demonstrating the metallicity enrichment of 
stars with planets with regard to c omparison "single" 
stars a nalysed in a homogeneous way l)Santos et al.ll200ll 
l2003h[ and references therein). This study provides pre- 
cise values of the effective temperatures, metallicities and 
gravity estimates, using a standard local thermodynam- 
ical equilibrium (LTE) analysis. These values have also 
been updated by using b etter recent oscillato r strengths 
(log g/) in the procedure (jSantos et al.ll2003a(l . From cal- 
ibrations of the width and surface of the CO RALIE cross- 
corre lation functions (CCF; described in ISantos et al.l 
|2002|) . we also derive the stellar projected rotational ve- 
locity vsini and [Fe/H] when not available from the spec- 
tral analysis. Like the majority of stars hosting planets, 
most of our candidates show a significant metal enrich- 
ment compared to the Sun, with the noticeable exception 
of HD 6434 (the second most deficient star hosting a planet 
known to date) and HD 111232. 

From the colour index, the measured Toff and the cor- 
responding bolometric correction, we estimate the star 



luminosities^ and we then interpolate the masses and 
ages in the grid of Geneva stellar evolutionary inodels 
with appropriate me tal abundances ijSchaller et alJll993 
ISchaerer et alJll993j) . For the brightest sta r in the sample 
we estimate, following ISantos et al.l l)2000|) . an activity in- 
dicator from the reemission flux in the Ca ii H absorption 
line. In some cases , this value is also available in the lit- 
erature (mainly in lHenrv et al.lll996|) . Such an indicator 
is then used to derive calibr ated estimates of the stellar 
rotat ional periods and ages l|Noves et al.l 1X984 lDonahu3 
Il993() . Tablengathers the photometric, astrometric, spec- 
troscopic information and inferred quantities available for 
our sample of stars hosting exoplanet candidates. 

The activity indicator is, however, not always available 
in the literature or cannot be estimated from our spec- 
tra because of the star faintness (HD 65216, HD 114386). 
Calcium reemission can then be visually checked on the co- 
addition of the CORALIE best-S/N spectra. For our sub- 
sample, this is shown in Fig. ^ for the A 3968.5 A CaiiH 
absorption line region. A prominent reemission feature is 
clearly visible for HD 114386 and HD 147513, the latter 
presenting, moreover, a strong activity indicator. However, 
these stars are only slowly rotating and no large influ- 
ence on the radial- velocity measurem ents is then expected 
ijSantos et al.ll200ot ISaar et alJll998l) . Hints of reemission 
are also observed for HD 65216, HD 83443, HD 142415, and 
HD 216770 but again usini is fairly low for these stars for 
which only a moderate jitter might be expected. 

Some stars in this sample deserve a few more com- 
ments: 



^ When not available from the sp ectroscopic study, a pho- 
tometric Tcff (following lFloweJll99a not given in Table0 is 
used to derive an indicative stax luminosity 



M. Mayor et al.: The CORALIE survey for southern extra-solar planets. XII 




3965 3970 3965 3970 3965 3970 
A(A) [Ca II H] 

Fig. 1. A 3968.5 A CaiiH absorption line region of the summed CORALIE good spectra for our candidate stars. The 
dotted lines indicate the exact position of the absorption hne centers. A clear reemission feature is visible for HD 114386 
and HD 147513, whereas only hints of reemission arc observed for HD 65216, HD 83443, HD 142415, and HD 216770. 
The spectra have been cleaned as much as possible from the pollution of the simultaneously recorded thorium-lamp 
reference spectra. 



- HD6434 (HIP 5054): The star is classified as sub- 
giant. Taking into account its low metallicity, the derived 
absolute magnitude, effective temperature and gravity es- 
timates also support the slightly evolved status of the star. 

- HD 19994 (HIP 14954, HR962, GJ128A): The star 
is known to have a close physical M dwarf companion 
(GJ 128 B) a few arcseconds away (- 100 AU; .Hal&.1994) . 
Planets in binary systems are important for our under- 
standing of planet formation because they seem to present 
different mass and orbital properties than planets orbit- 



ing si n gle stars JZucker fc Mazehl2002HEggenberger et alJ 
l2003t lUdrv et alJ l2003d|) . The rotational velocity of 
HD 19994 is fairly large (w sinz = 8.1 kms""'^) thus, even if 
the star is not clearly active, some small radial- velocity jit- 
ter might be expected. Contamination from the secondary 
could also be a concern although the 6.4 magnitude dif- 
ference between the stars makes it very unlikely. 

- HD 82943 (HIP 47007): A fair amount of the easy- 
to-burn elem ent ^Li ha s been detected in the spectra of 
this star Ijisraelian et al...200L .2003.) suggesting a planet 
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engulfment in the late stages of the system formation, af- 
ter the star has reached the Main Sequence. The not too 
deep convective zone of this GO dwarf allows then for the 
survival of this element. 

- HD 111232 (HIP 62534): The high velocity 
(y = 104.4kms-i) and low metallicity ([Fe/H] = -0.36) 
of the star indicate that it probably belongs to the thick 
disk population. HD 111232 is a proposed binary in the 
Hipparcos catalogue It is one of the so-called 
problem stars for which an astrometric acceleration 
solution is provided (flagged "G" in the H59 field of the 
main catalogue). However, no companion with Ay < 3.0 
was found clos e to the star (within 1.08") by speckle 
interferometry ijMason et al ]|l998). 

- HD 121504 (HIP 68162): A visual companion 

i CPD —55:5793) is obser ved at a separation of 34.2" 
Dommaneet fc Nvg|ll99l . It is an A2 star of magnitude 
y = 9.17 with different proper motions than HD 121504. 
The pair is thus not physical. No compani on was found 
close to the star by speckle interferometry ( Mason et al.l 
11998^ . 

- HD 142415 (HIP 78169): A Rosat- All-Sky-Survey 
X-ray source (IRXS J155740.7-601154) is observed close 
to the star, at ~5". No bright companion was found 
close to the star by speckle interferometry measurements 
()Ma,son et a,lJll99Sl) . 

- HD 147513 (HIP 80337, HR6094, GJ 620.1 A): 
The star was proposed to b e a barium dwarf by 
IPorto de Mello fc da Silval ^ll997^ . A common proper mo- 
tion white dwarf (WD; y = 10), at a projected distance 
of 5360 AU, could support the explanation of the ori- 
gin of the barium feature by the process of mass trans- 
fer in a binary system, in which the secondary com- 
ponent accreted m atter from the evoly ed primary, now 
the WD (see e.g. I.Torissen et all Il998l and references 
therein). To account for the ob served large separation, 
IPorto de Mello fc da Silval l|l997^ invoke a possible ejec- 
tion of the WD from an originally quadruple system also 
including HD 147513. In such a case, this would be the 
first known case of a planet orbiting a solar-type star in 
a binary system with the stellar companion in the last 
stage of its life. The influence on the planet evolution of 
the mass transfer between two stars is still poorly stud- 
ied. The mentioned authors also include the star in the 
0.3 Gyr old Ursa Major kinematical group. Although some 
doubts (but no r ejection) were cast on the membership 
()King et a,].ll20n3l) . it is worth noticing that the activity 
level, activity-derived age and metallicity (Tabled corre- 
spond well to the Ursa Major group. Finally, no bright 
companion was found closer -in to the sta.r by s peckle in- 
terferometry measurements l)Mason et al.lll99^ . 

3. Radial-velocity data and orbital solutions 

The radial-velocity measurements presented in this pa- 
per were obtained with the CORALIE echelle spectro- 
graph mounted on the 1.2-m Euler Swiss telescope at 
La Silla Observatory (ESQ, Chile). CORALIE is a simi- 



lar but improved version of the ELODIE spectrograph at 
the Haute-Provence Observatory (CNRS, France). Details 
on the instrument design, reduction procedures as well as 
radial-velocity estimates based on simultaneous thorium- 
lamp meas urements and cross-co rrelation technique can 
be found in lBaranne et all lll996l) . 

The typical precision obtained with CORALIE is 
~3ms~^ for bright stars (Oueloz et al. 2001b}. However, 
due to the small s ize of the telescope, our sample stars 
ljUdrv et al.ll200o'ah are mainly photon-noise limited. For 
the derivation of the orbital solution an instrumental Bi- 
tot of 3 m s~^ is quadratically added to photon noise. Then 
we only take into account the spectra with good signal- 
to-noise, typically corresponding to radial-velocity uncer- 
tainties below 10 to 15ms~^, depending on the star. Also, 
to improve the Doppler information extraction from the 
spectra, a new weighted cr oss-correlation scheme was de- 
veloped ijPepe et al.ll2002aj) . It is applied for the velocity 
estimates of most of the stars of this paper^. 

3.1. The single-planet systems 

To simplify the presentation of our results for the single- 
planet systems, the standard orbital parameters derived 
from the best fitted 1-Keplerian solution to the data are 
gathered in Table|21 Useful inferred planetary parame- 
ters and further information on the number of measure- 
ments used, their time coverage (Span) and the resid- 
uals (weighted r.m.s.) around the solutions (cr(0 — C)) 
are given as well. Finally, we also indicate the refer- 
ences for the detection announcements of these candi- 
dates^ {Det.Ref.). The corresponding phase-folded radial- 
velocity curves are displayed in Figs. [3 to 0] (top panel in 
each diagram). The residuals drawn as a function of the 
Julian date are displayed in the figures as well (bottom 
panels). 

Simultaneous independent d iscoveries have been m ade 
for 3 of our ca ndidates: by iFischer et 

ID l|200lh for 

HD 92788 and by I.Tones et all l|2002t) for HD 196050 and 



The CORALIE numerical mask used in the cross correla- 
tion for the velocity estimate was built from a KO-dwarf spec- 
trum. For much earlier spectral-type stars (typically < Gl), 
the mismatch between the stellar spectrum and the template 
is enhanced by the weighted scheme and no improvement is 
obtained for the radial-velocity measurements. The usual cross 
correlation scheme is then used in such cases (Table|5J 

^ At the lAU Symp. 202 in M anchester, we announc ed 6 
CORALI E candid ates (iMavor et al. 2000; Oueloz ct al. 2()0ol : 
fUdrv et all l2000bl 5 of them being described in th is pa- 
per) and 1 further ELODIE candidate (ISivan et alJl2000h . The 
proceedings of the lAU Symp. 202 have not yet appeared. 
Our 8 contributions to this conference describing the planet 
and brown-dwarf CORALIE and ELODIE new detections, the 
metallicity of stars hosting planets, the effect of stellar ac- 
tivity on radial-velocity measurements and a presentation of 
the new HARPS spectrograph (now running on the ESO 3.6- 
m at La Silla) are accessible from our Exoplanets web page: 
|http://obswww.unige.ch/Exoplanets/publications.html, 
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Table 2. CORALIE best Keplerian orbital solutions as well as inferred planetary parameters for the 1-planet systems. 
Span is the time interval in days between the first and last measurements. (t(0 — C) is the weighted r.m.s. of the 
residuals around the derived solutions, mask is the template used in the cross-correlation scheme for the radial-velocity 
estimate (see Sect. 3). Det.Ref. is the reference to the planet detection announcement (conference or press release). 
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-2.278 ±0.004 31.153 ±0.002 




Ul 


[deg] 


147 ±12 


63 ±22 


281 ± 10 




K 


[ms-i] 


55.0 ±6.2 


34.6 ±5.7 


30.9 ±1.9 




ai sin i 


[10-3 AU] 


6.367 


3.815 




0.313 




fim) 


[10-9 Ms] 


19.74 


4.694 




0.291 




m2 sin 2 


[Mj,p] 


3.02 


1.82 




0.65 




a 


[AU] 


2.43 


2.32 




0.46 




^ mcas 




31 


21 




16 




Span 


[days] 
[ms-i] 


1364 


1405 




827 




a{0 - C) 


7.2 


7.2 




7.8 




mask 




weighted KO 


weighted KO wei; 


ghted KO 





Independently discovered bv lFischer et al .1 ll200J) f" Independently discovered bv lJones et al.l \200± 
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HD 216437. Their orbital solutions are very similar to ours 
(TableEl. 

For most of the systems, the residuals measured 
around the derived solutions are compatible with the typ- 
ical individual radial- velocity uncertainties. Some of them 
are, however, slightly larger (~10ms~^). In these cases, 
a slight level of activity can be invoked to explain the 
additional noise. HD 114386 is a good example with its 
clear activity level testified by the visible reemission in 
the CaiiH absorption line (Fig.^. Another illustration 
is given by the similar early G dwarfs HD 121504 and 
HD 142415. They present light rotation and an activity 
level that can explain the somewhat large measured resid- 
uals. It is interesting to note here the similar values of 
vsini, logR'jjj^ and a{0 — C) for the two candidates. 

In order to emphasize a relation between the residuals 
to the Keplerian solutions and stellar activity, we com- 
pared in a systematic way for all the candidates the ob- 
tained residuals with the shape of the spectral lines, esti- 
mated by the bis ector inverse slope of our cross-correlation 
functions CBIS: lOueloz et all l2001a|) . Unfortunately, at 
this level of velocity variation (<10ms~^) and at the 
S /N of the CORALIE spectra, no definitive conclusion can 
be drawn. The same non-conclusive result was also ob- 
tained for HD 73256 which did not show any clear rela- 
tion between the BIS parameter and the residuals around 
the 2.5-d solution, although the activity- induced radial- 
velocity jitter was undubitably emphasized by the simul- 
taneous varia tions of the stellar photometric signal and 
the residuals l|Udrv et al.ll2003a| . We also searched in a 
systematic way, through Fourier analysis, for additional 
periodicities in the residuals around the derived solutions. 
Nothing significant was found for these candidates. Some 
of them deserve, however, further comments (see the fol- 
lowing subsections). 

Amongst the 13 single-planet systems presented in 
this subsection, 10 (HD 19994, HD 65216, HD 92788, 
HD 111232, HD 114386, HD 142415, HD 147513, 
HD 196050, HD 216437, HD 216770) are on interme- 
diate or long-period orbits (100 < P < 1350 days) with 
medium eccentricities (0.2<e<0.5). Such properties 
are shared by the bulk of the presently known exoplan- 
ets^. The remaining 3 shorter-period planets (HD6434, 
HD 83443, HD 121504) present, on the other hand, a 
lower eccentricity. The eccentricity-period trend is also 
observed for stellar binaries. This similitude between 
the two populations was often brought up as a question 
for different planet and binary for mations. Nevertheless, 
clear differences exist (see e.g. iMavor fc Udrvl l2000t 
[Halbwachs et al. 2003, for a more complete discussion). 
Also the mass-period relat ion of exoplanets , checke d for 
statistical significan ce bv [ Zucker fc Mazehl (jifl^^) ^^id 
further discussed in lUdrv et aL rT2003d), is apparent in 

^ Long-period candidates on low-eccentricity orbits (closer 
to the giant planets in our solar systems) are however more 
and more often discovered as the time coverage of the radial- 
velocity surveys increases 



our subsample; the two shortest-period systems host (by 
far) the lightest planets, with minimum masses below 
0.4Mjup, whereas the longest periods (P> 1000d) are 
associated with planet minimum masses of 1.82, 3.02 and 
6.8Mj^p. 

Comments on HD6434 

Since the discovery of the planet, we have doubled the 
number of measurements available for HD 6434, gathering 
a total of 130 good spectra over more than 1500 days. The 
star relative faintness and low metallicity lead to a typi- 
cal radial-velocity photon-noise uncertainty on individual 
measurements of only ~8ms~^. However, the large num- 
ber of measurements allows us to derive precise Keplerian 
orbital elements for the system (Table|2l. The very low 
planet minimum mass inferred from the orbital solution 
(rn2 sin i = 0.39 Mjup) gives us a first example of a very 
light planet orbiting a deficient star. 

Although the star is found to be neither active nor 
rapidly rotating, a concern is brought by the not so dif- 
ferent values of the rotational period (18.6 day) estimated 
from the activity indicator and the orbital period (22 day). 
The uncertainty on the orbital period is very small, but the 
calibrated P^ot carries an intrinsic uncertainty. The a pos- 
teriori verification that the shape of the spectral lines, esti- 
mated by the bi sector inverse slope o f our cross-correlation 
functions fBlS: lOueloz et al.ll2001a[> . does not vary with 
the radial velocity provides an indication against activity 
to be the source of the observed radial- velocity variation. 
The BIS itself is even constant with a r.m.s. of 8.6 ms~^, at 
the photon-noise level. Moreover, the radial-velocity vari- 
ation is stable over more than 68 c ycles. However, rough 
simulations bv lSantos et~al] l)2003(Jl raised the possibility 
of radial-velocity variations without noticeable change in 
the BIS for slowly rotating stars. In such cases, only simul- 
taneous velocity and photometric measurements can trace 
the intrinsic ori gin of the radial-velo city variations, as e.g. 
for HD 19 2263 llSantos et al.ll2003(J) or for the HD 73256 
residuals l|Udrv et^dll2003^ . No indication of a 22-d pe- 
riodicity is present in the Hipparcos photometric data for 
HD 6434 and the Geneva photometry finds the star stable 
at a 3 mmag level. 

Comments on HD 19994 

The residuals around the Keplerian solution (8.1 ms^^) 
are slightly larger than individual photon-noise errors 
(median at 6ms~^). They seem mainly due to veloci- 
ties around JD = 2 451600 presenting a large dispersion. 
Taking into account the fairly large rotation of the star 
and its binarity status, this could possibly come from an 
enhanced stellar activity level at that moment, although 
the BIS parameter does not show anything special at that 
time. Nothing else particular is visible in the residuals. 
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Fig. 4. Phase-folded radial-velocity measurements ob- 
tained with CORALIE for HD 83443, superimposed to the 
best Keplerian planetary solution (top). The residuals as 
a function of Julian Date are displayed in the lower panel. 



Comments on HD 83443 

HD 83443 was first announced at the Manchester lAU 
Symp. 202 to host a resonant 2-plai iet system with pe- 
riods Pi=2.985d and P2=29.85d ijMavor et al. 20001. 
At that time the available 93 high-S/N CORALIE spec- 
tra allowed us to clearly find the 2 periods of variation. 
The two periodic signals were highly significant. In par- 
ticular, the false alarm probability for the 30-d period 
in the data was very low (<9 • 10"''), close to the 4cr 
detection limit (Fig.[3 top). Moreover, these periodici- 
ties were not present in the velocities of contemporary 
observed constan t stars with similar photometric charac- 
teristics llUdrv et al...2002 '). This rules out instrumental 
effects as the source of the observed variations. A new re- 
duction of the data using the weighted cross-correlation 
scheme also confirms the results presented in Manche ster. 

However, about 2 years later. iButler et all lj2002l) try- 
ing to confirm the 2-planet system did not detect the 30-d 
second-planet signal in their Keck-|-AAT radial-velocity 
data. We then also verified that it had disapp eared in our 
more recent data as well (jUdrv et al.ll20o3) . The corre- 
sponding peak in the Fourier transform of the radial ve- 
locities is no longer present (Fig.El bottom). 

The origin of this transient signal is not clear yet. An 
appealing possibility is to attribute the effect to activity as 
the corresponding period is compatible with the activity- 
related stellar rotation estimate (within the uncertainties 
related to the calibration). As mentioned above, we tried 
for this star as well to see a relation between BIS and the 
residuals around the 2.98-d Keplerian solution but without 
success at the precision of our CORALIE data. No 30-d 
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Fig. 5. Fourier transform of the residuals around the 2.98- 
d orbital solution of HD 83443 b Top: for radial- velocities 
in the 2 451 500 < JD < 2 451 702 period, when the 30-d pe- 
riod signal was seen in the data, and Bottom: for later 
velocities. The false alarm probability of the 30-d period 
in the early measurements is only 9 • 10"** (close to the 4 a 
detection level) , whereas the signal has disappeared in the 
later data. 



periodicity is found in the BIS data either. However, if 
existing, such a relation between residuals and line shape 
should easily come out of the HARPS data that are already 
providing a preliminary commissioning or bital solution for 
HD 83443 at a - 1.5 ms^^ precision level l|Mavoill20ol . 

Looking at the orbital solution given in Tabled we 
can note that the ecccentricity derived for the planet is 
not significantly different from 0. 

With its 2.98-day period, HD 83443 b was a good candi- 
date for photometric transit search. The photometric uvby 
observations were obtained with the Stromgren Automatic 
Telescope (SAT) at ESO LaSilla, Chile (Olsen et al. in 
prep). Unfortunately, no transit was detected. 
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Table 3. CORALIE simulataneously derived 2-Keplerian orbital solutions as well as inferred planetary parameters 
for the multi-planet systems HD 82943 and HD 169830. Parameter definitions are the same as in Table|5] Note that for 
HD 82943, the e and oj elements are not well defined in the sense that there exist other solutions with very different e 
and Lo estimates but with equivalent <t(0 — C) values (see text). 



Parameter 


HD 82943 c 


HD 82943 b 


HD 169830 b HD 169830 c 


npf ftp f 




(HIP 47007) 
ESO (May 2000); ESO (April 2001) 


Naef et al. ( 


(HIP 90485) 
2001): Udrv et al. ( 2003b) 




-J- ^ 

[daysj 


o 1 n /I _i_ n o 


/I Q C 1 _l_ 1 A 


225.62 ±0.22 2102 ±264 


T 

J. 


[JU-i ^UU UUUJ 


52284 ± 1 


51758 ±13 


51923 ±1 


52516 ±25 


G 




0.38 ±0.01 


0.18 ±0.04 


0.31 ±0.01 0.33 ±0.02 


V 


[kms"-^] 


8 


144 ±0.001 




-17.209 ±0.006 




[deg] 


124 ±3 


237 ±13 


148 ±2 


252 ±8 


K 


[ms-^] 


61.5 ±1.7 


45.8 ±1.0 


80.7 ±0.9 


54.3 ±3.6 


ai sin i 


[10-^ AU] 


1.145 


1.801 


1.591 


9.898 




[10-3 Mq] 


4.156 


4.118 


10.56 


29.27 


1712 sini 


[Mjup] 


1.85 


1.84 


2.88 


4.04 


a 


[AU] 


0.75 


1.18 


0.81 


3.60 


^ mcas 






142 




112 


Span 


[days] 




1593 




1506 


a{0 - C) 


[ms-i] 




6.8 




8.9 


mask 




weighted KQ 




weighted KQ 



Comments on HD 121504 

We mentioned above that the somewhat large residuals 
obtained for this system may probably be related to the 
stellar activity-induced jitter (logi?^^- = —4.57). We also 
could see in the temporal distribution of the residuals some 
indications of an additional radial-velocity drift. If real, 
this drift is very small. A combined Keplerian ± linear 
drift model yields a drift value of '^3ms-^y-^, with- 
out changing noticeably the planetary orbital parameters. 
Furthermore, the combined fit does not improve much the 
quality of the solution: the a{0 — C) decreases only from 
11.6 to 11.2ms^^. The drift is thus not considered as sig- 
nificant and not included in the solution given in Table|5] 
Future measurements will confirm or rule out this poten- 
tial drift. 

Comments on HD 142415 

The Keplerian solution derived for HD 142415 yields a pe- 
riod of 386.3 ±1.6 days, not too far from 1 year. As we 
have points covering the maximum and close (on both 
sides) to the minimum of the phased radial- velocity curve 
over 4 full cycles, this proximity to an annual variation 
is not so much of a concern for the reality of the plane- 
tary signal. The trouble resides rather in the difficulty of 
covering the whole phase interval. As a result the mini- 
mum of the curve is highly undersampled and the best 
fitted orbital solution tends to be very eccentric, increas- 
ing the radial-velocity semi-amplitude in an artificial way. 
This effect is even worsened by the stellar activity induced 
jitter that makes the orbital elements more difficult to de- 
termine. We have thus arbitrarily fixed the eccentricity of 
the solution to e = 0.5, a plausible value when looking at 
the sequence of the radial-velocity measurements. Note, 



however, that solutions with e between 0.2 and 0.8 would 
still be acceptable. 

Comments on HD 196050 

The phase coverage of our 31 CORALIE radial veloci- 
ties is not good enough to obtain a constrained value of 
the planet eccentricity. We thus fixed this eccentricity to 
e = 0.3, a value minimizing the residuals around the de- 
rived solution. 

3.2. The 2-planet systems 

3.2.1. The 2-planet resonant system around HD 82943 

The star hosts a very interesting 2-planet system. Between 
JD = 2 451184 and JD = 2 452 777, we have gathered 142 
good CORALIE spectra of the star that allow us to 
derive a 2-Keplerian approximate solution with periods 
Pb = 435.1 d and Pc = 219.4 d. The whole set of orbital pa- 
rameters are given in Table|31and the solution and veloci- 
ties are displayed in Fig.|Hl (left). The particular values of 
the periods made us miss the short-period minimum of the 
curve when the star was behind the Sun. The long-perio d 
planet (HD 82943 b) was thus announced first ||ESQ||2000|) . 
about 1 year before the detection of the secondary planet 
signature (HD 82943 c; EgO 20Q1). 

The Pi /Pj = 2/1 resonant systems are very impor- 
tant because, when the planet orbital separations are not 
too large, planet-planet gravitational interactions become 
non-negligible during planet "close" encounters, and will 
noticeably influence the system evolution on a timescale 
of the order of a few times the long period. The radial- 
velocity variations of the central star will then differ sub- 
stantially from velocity variations derived assuming the 
planets are executing independent Keplerian motions. We 
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observe a temporal variation of the instantaneous orbital 
elements. In the most favourable cases, the orbital-plane 
inclinations, not otherwise known from the radial- velocity 
technique, can be determined since the amplitude of the 
planet-planet interaction directly scales with their true 
masses. 

In the case of multiple planets, only approximate ana- 
lytic solutions of the gravitational equations of motion ex- 
ist, and one must resort to numerical integrations to model 
the data. Several studies have b een conducted in this 
directio n for the G1 876 system ('Laughlin fc Chamber j 
[200 L ..Rivera fc Lissauer 2001: Nau cnbcrg 200 similar 
to HD 82943 in the sense that it also hosts two 2/1 res- 
onant planets at fairly small separations. The results of 
the Newtonian modeling of the G1876 system have vali- 
dated the method, improving notably the determination 
of the planetary orbital elements, but the time coverage 
of the measurements is still too small for the method to 
provide strong constraints on the plane inclinations. The 
valley of the acceptable solutions is still very shallow, al- 
though including the cor rect answer pro vided by HST 
astrometric observations ((Benedict et al.l 12002'). Further 
radial-velocity measurements will undoubtedly improve 
the situation. 

Our correct modeling of HD 82943, taking into account 
the planet-planet interactions, is under study and will be 
presented in a forthcoming paper with a more detailed de- 
scription of the system behaviour (Correia et al. in prep). 

Our present approximate solution for the system 
(TableO yields residuals at the level of the photon noise of 
the radial- velocity measurements. The inferred two plane- 
tary masses are very similar. Although in a 2/1 resonance, 
the two orbits do not seem to be aligned. However, we have 
to warn the reader not to take too litterally the results de- 
scribed here because: 

- As mentioned above, due to the planet-planet inter- 
actions, the orbital elements are time dependent; 

- Due to the non-optimum phase coverage, the planet 
eccentricities and ws are badly (if at all) constrained by 
the data. The small uncertainties given in the table only 
relate to the given local solution in parameter space. There 
are, however, other local solutions with almost equivalent 

minimum values. For instance, a solution with the very 
different values of the eccentricities Cc = 0.4 and Cf, ~ only 
increases the residuals from 6.8 to 6.9 ms~^. In this case, 
as ef, ~ 0, LUb is completely non-determined and it is then 
possible to find for these eccentricities an aligned configu- 
ration with Wf, = Wc = 110 deg and the same level of resid- 
uals (6.9 ms~^), leaving furthermore Pi and Ki almost 
unchanged. 

We will now continue to follow closely this system, 
accumulating measurements to improve in the future the 
derived solution. 



3.2.2. The HD 169830 hierarchical 2-planet system 

A 230-d period plaii e t orb iting HD 169830 was first de- 
scribed in iNaef et alJ ll200lh. The orbital solution was de- 
rived from the 35 CORALIE spectra available at that time. 
After the second maximum of the radial- velocity curve, we 
noticed an additional trend in the data pushing us to fol- 
low the star more closely. We have now gathered 112 good 
spectra that allow us to simultaneously derive a complete 
2-Keplerian orbital solution for the system (TableOl and 
Fig.El right). 

Unlike HD 82943, this system is not resonant but more 
hierarchically structured. The separation between the 2 
planets stays always fairly large, the 2-Keplerian model 
is then supposed to provide a good approximation of the 
system evolution on a fairly long time. However, the time 
span of our velocity measurements barely covers the long 
period variation. The corresponding planetary orbit is 
thus not completely constrained and a large uncertainty 
is still observed for the long period. 

This 2-planet system is the first to be an nounced after 
the proposition by Mazch fc Ziickeil ( 2003bl) of a possible 
correlation between mass ratio and period ratio for adja- 
cent planets in multi-planet systems. It is interesting to 
note t hat the new system agree s with the proposed corre- 
lation 1 Mazeh fc Zuckeif 



4. Summary 

We have described in this paper 16 still unpublished exo- 
planet candidates discovered with the CORALIE echelle 
spectrograph mounted on the 1.2-m Euler Swiss telescope 
at La Silla Observatory. Amongst these new candidates: 

- Ten are typical extrasolar Jupiter-like planets 
on intermediate- or long-period (100 < P < 1350 d) and 
fairly eccentric (0.2 < e < 0.5) orbits (HD 19994, HD 65216, 
HD 92788, HD 111232, HD 114386, HD 142415, HD 147513, 
HD 196050, HD 216437, HD 216770). They resemble the 
bulk of extra-solar planets found to date. 

- Two of these planets (HD 19994, HD 147513) are 
orbiting one component of a multiple-star system. 
Such planets seem to present different orbital and 
mass characterist i cs th a n the other singJe - star p lanets 
((Ziicker fc Mazehl l2002t lEggenberger et all l2003ll . The 
companion to HD 147513 is even a white dwarf, the evo- 
lution to which has probably also influenced the planet 
evolution through mass transfer between the two stars. 

- Three candidates are shorter-period planets 
(HD6434, HD 121504, HD 83443) with lower eccentricities 
(the latter being a hot Jupiter). 

- More interesting cases are given by the multiple- 
planet systems HD 82943 and HD 169830. HD 82943 is a 
resonant Pt / Pc = 2/1 system in which planet-planet inter- 
actions are influencing the system evolution. HD 169830 
is non-resonant and more hierarchically structured, and 
therefore less affected by this kind of interaction. 

From a more global point of view, our candidates fol- 
low the period-eccentricity and period-mass trends ob- 
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served for the whole sample of known extra-solar planets. 
They follow as well the trend for stars hosting planets to 
be more metal rich than normal stars of the solar neigh- 
bourhood llSantos et al.l200l[|2003btlGonzalez et al.l200ll: 
iLaws et al.ll2003l) . Only 3 amongst the 15 stars are metal 
deficient with regards to the Sun whereas almost all the 
others present high [Fe/H] values. 

We emphasize the difficulty encountered to fully con- 
strain multi-planet systems. Such a task, involving many 
free parameters, requires a good phase coverage and a 
fair number of measurements, even for the simplest cases. 
As a consequence, studies on multi-planet system stabil- 
ity should not rely too closely on the given orbital pa- 
rameters. The published solutions will probably change 
in the future (some will notably change) as more mea- 
surements become available. A substantial advance in this 
domain will be brought by the the new HARPS spectro- 

j raph mounted oii the ESO 3.6-m telescope at La Silla 
Pepe et alJ f2002bj) available since October 2003. With 
the very high precision achieved for radial-velocity mea- 
surements and the quality of the spectra, HARPS is now 
providing us with an unequalled tool to characterize multi- 
planet systems and/or disentangle activity-induced jitter 
from orbital radial- velocity variations. 
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Fig. 2. Phase-folded radial-velocity measurements obtained with CORALIE for HD6434, HD 19994, HD 65216, 
HD 92788, HD 111232 and HD 114386, superimposed on the best Keplerian planetary solution (top panel in each 
diagram). The residuals as a function of Julian Date are displayed in the lower panels. 
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Fig. 3. Phase-folded radial-velocity measurements obtained with CORALIE for HD 121504, HD 142415, HD 147513, 
HD 196050, HD 216437 and HD 216437, superimposed on the best Keplerian planetary solution (top panel in each 
diagram). The residuals as a function of Julian Date are displayed in the lower panels. 
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Fig. 6. Top: Temporal radial-velocity measurements obtained with CORALIE for the 2-planet systems HD 82943 (left) 
and HD 169830 (right), superimposed to the best simultaneously-derived 2-Keplerian planetary solutions. Bottom: 
Residuals as a function of Julian Date. 



